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It has been clear for the past several years that new physics in the quark sector can only appear, 
in low energy observables, as a perturbation. Therefore precise theoretical predictions and precise 
experimental measurements have become mandatory. Here we draw attention to the significant 
advances that have been made in lattice QCD simulations in recent years in K ^ tttt, in the 
long-distance contribution to indirect CP violation in the Kaon system (e) and in rare A-decays. 
Thus, in conjunction with experiments, the construction of a unitarity triangle purely from Kaon 
physics should soon become feasible. We want to emphasize that in our approach to the A-unitarity 
triangle, the ability of lattice QCD methods to systematically improve the calculation of the direct 
CP-violation parameter (s') plays a pivotal role. Along with the B-unitarity triangle, this should 
could allow, depending on the pattern of new physics, for more stringent tests of the Standard Model 
and tighter constraints on new physics. 


INTRODUCTION AND MOTIVATION 

For the past decade or more, from a variety of low 
energy precision experiments, such as those from B- 
factories and LHCb, it is becoming clear that effects of 
new physics are likely to just show up as perturbations. 
The current tests on the unitarity triangle show very 
good agreement with the Standard Model (SM) CKM- 
paradigm mm to the level of about 20%. Impressive 
as this is, it is also important to emphasize that 10-15% 
effects of new physics (which at present can not be ruled 
out) are quite big. In fact it is useful to recall that CP 
violation was first discovered in the decays Kp ^ tttt 
at the level of 10“^ [3]; therefore we should not be sur¬ 
prised if similar precise measurements become necessary 
to discover new phenomena. Thus the need for improved 
precision in experiments as well as in theory should be 
clear. 

In this work we highlight recent advances made in the 
lattice computation of A-decays that had been serious 
challenges for a very long time. As will be explained, 
in the course of tackling the calculation of the crucially 
important direct CP violation parameter in A ^ tttt 
decays, e'/s, by using the Lellouch-Liischer method in¬ 
volving finite-volume correlation functions, the RBC and 
UKQCD collaborations managed to develop an interest¬ 
ing method to tackle matrix elements of non-local 4- 
quark operators which are relevant to quantify “Long- 
Distance” (LD) contributions to a variety of matrix ele¬ 
ments of physical interest. Explicit examples under cur¬ 
rent active study are the neutral Kaon mass difference, 
the related long-distance part of the indirect CP-violation 
parameter 5 and the branching ratio (BR) of the rare de¬ 
cay A+ ^ TT+Z/P. 

For e'/s, e, and Atuk rather precise experimental mea¬ 
surements already exist [4]. The initial measurement of 
BR(A+ ^ TT+z^P) = (l.TStl’.Js) 10“^^ was accom¬ 
plished some years ago at BNL in experiment E-787 and 
E-949 [5l|6]. At CERN the NA62 experiment is expected 


to significantly improve the determination of this branch¬ 
ing ratio [31] in the next few years. 

An important point to consider is also that lattice 
methods for tackling non-perturbative effects are largely 
systematically improvable. This means once the physical 
quantity becomes amenable to lattice methods, accuracy 
with computer capability and time is essentially guaran¬ 
teed. A few examples are the Kaon A-parameter Bk 
(that is the ratio of the complete A-A matrix element 
of the leading SM operator to its estimate in the naive 
factorization approximation) which in full QCD with chi¬ 
ral fermions got evaluated with about 7-8% total error 
around 2007. In the next ^ 5 years, many different 
collaborations have attacked it with different discretiza¬ 
tions and the world average now has an error around 1%. 
Quark masses, A ^ tt^z/, A ^ tttt in the 1 = 2 channel, 
/^, A-mixings, and others have followed a similar path 
(see, e.g., Ref. m for a review of the present status of 
flavor physics lattice calculations). Calculations of e'je 
and the long-distance effects in A-A mixings and in rare 
A decays are expected to progress in analogous ways. 

Bearing in mind these exciting developments in theory 
and in experiments, it seems timely to ask if we can now 
construct a unitarity triangle based primarily from input 
from A-physics, which has been often talked about m 
[16], and is the subject of this work. 

By way of motivation, let us recall that, in general, nat¬ 
uralness arguments suggest that beyond the SM (BSM) 
scenarios are unlikely to be flavor blind. Indeed, just as in 
the SM, as weak interactions are “switched-on” the gauge 
eigenstates are no longer aligned with the mass eigen¬ 
states and the connection between the two is monitored 
by the CKM-matrix. Warped extra dimensions provide 
a very interesting example, as they provide a geometric 
understanding of flavors. Many aspects of flavor physics 
can be readily understood through localization of differ¬ 
ent flavors at different locations in the extra-dimension 
yielding a non-universal, strongly mass-dependent effect 
on flavors accompanied by many 0(1) BSM-phases pTL 
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US]. This is a very illustrative example, and in many if 
not most BSMs a similar situation arises. 

Another important consideration by way of motivation 
is that the observables in the iC-unitarty triangle are 
widely believed to be very sensitive to deviations from 
minimum flavor violation (MFV); this is in contrast to 
the many observables used in the standard 5-unitarity 
triangle. Thus using the B-UT to extract SM-CKM pa¬ 
rameters and then the iC-observables for new physics 
searches can be a very effective approach. 

The ability of these new lattice methods to quantify 
non-perturbative effects in K-ir physics calls for a re¬ 
examination of Kaon experiments. For one thing given 
that in a few years the theory errors on s' are likely to 
come down to around 11 % of the current experimental 
central value with the real possibility of further improve¬ 
ments down the road, improved experimental determina¬ 
tions may be called for given the current experimental 
error is around 15%. 

Moreover, although measurement of s' is very chal¬ 
lenging, measuring the rate for Kl ^ is even more 

challenging. Traditionally, this very difficult purely CP- 
violating rare mode m has been a pristine SM predic¬ 
tion as it is clearly short-distance dominated because of 
the large top quark mass. However, as lattice meth¬ 
ods make sufficient progress in s' ^ the theoretical advan¬ 
tage of Kl will likely diminish. Note, however, 

that contributions to these two processes are indepen¬ 
dent in many models making both invaluable tools in 
searches for new physics effects. Moreover, it must be 
stressed that progress towards a precise determination of 
BR(i^+ ^ 7r~^uu) is guaranteed by existing and planned 
experiments and theoretical progress in lattice QCD cal¬ 
culations of such processes (This will be achieved by re¬ 
duced errors on | Vcb \ from improved determinations of the 
B ^ D form factor (See for instance Ref. [ 21 ] where the 
dependence of this form factor is calculated) and by 
a direct calculation of long distance charm contributions 
to BR(i^+ ^ TT+z/p) [ 22 |). Needless to say, improve¬ 
ments on the measurement of s' js will require a major 
experimental effort. 

Another very interesting decay mode that, in principle, 
can provide a useful test of CP-violation is ^ 7 r^e+e“ 
(see, e.^.. Ref. [23]) which in fact allows for a detailed 
study via a Dalitz plot. The major theoretical distinc¬ 
tion with Ki^ is that the charged lepton mode 

receives a non-negligible contribution to the BR from CP- 
conserving 2 -photon intermediate state, as opposed to 
Kl, 'K^vv which, to an excellent approximation in the 
SM is purely CP-violating. 

If theory could reliably and precisely predict the CP- 
conserving contribution to Kj, 7 r^e+e“ which is dom¬ 
inated by LD effects, this mode could become extremely 
significant as from an experimental perspective the e+e“ 
final state is seemingly a lot easier to address. Unfor¬ 
tunately there is a daunting experimental background. 


pointed out in [24] from Kl ^ 770 +e“ which will need 
to be overcome [25] . 

The measurements that we include in the Kaon unitar- 
ity triangle fit (KUT) are 0 , s'/s, BR(i^+ ^ and 

\V^^\ (from inclusive and exclusive semileptonic b c£u 
decays). The standard unitarity triangle fit (SUT) de¬ 
pends on \ Vcb\ only via the ratio \ Vub/Vcb\ and, indirectly, 
via the rare decay B rz^^-. Presently \Vcb\ provides a 
subdominant contribution to the uncertainties on these 
two quantities (the former is controlled by theoretical 
and experimental errors on 6 ^ u£u decays, the latter is 
dominated by experimental errors). Hence the SUT and 
KUT fits are essentially independent and any tension be¬ 
tween them is an indication of BSM physics. Finally we 
note that Am^, while being very sensitive to the chiral¬ 
ity of new physics [26], does not provide any constraint 
on the p — fj plane and is therefore not included in the 
KUT fit. 

The structure of this manuscript is as follows: We first 
briefly review new lattice methodology that enabled some 
of the recent advances in lattice Kaon physics. We then 
specifically summarize K ttuP decays, s'/s, and 0 . 
We Anally present our findings for current and potential 
future Kaon unitarity triangle fits. 


LATTICE METHODOLOGY 


Relations between finite-volume Euclidean correlation 
functions, which are accessible to lattice QCD, and 
infinite-volume matrix elements form the foundation of 
many of the lattice efforts mentioned in this work. The 
simplest case of one particle decaying into a single two- 
particle final-state allows for the extraction of K ^ tttt 
matrix elements [23 EH). Extensions for matrix elements 
beyond this limit have recently received a lot of atten¬ 
tion, see Ref. [29] for a summary. Eor the current work, 
the development of methodology for non-local (bi-local) 
matrix elements is of particular interest [3QH32] . The Kl 
- Ks mass difference computation can serve as a nice in¬ 
troduction to the general methodology m- The desired 
infinite-volume quantity 


^ mx — En 


can be related to its finite-volume counter-part 
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with controlled finite-volume errors, where it is assumed 
that a single tt-tt intermediate state no is degenerate with 
K^ and K . The remaining task is to extract Am^ 
from Euclidean space correlation functions which is usu- 




3 


ally formulated in terms of four-point functions 

A=^(K\tf) p dh p dhHw{t2)Hw{ti)K°{ti)). 

(3) 

This amounts to the summation of operator insertions 
within a fiducial volume bounded by time coordinates 
ta and ti) with source and sink operators at U and tf 
satisfying tf ta < h ti. Inserting a full set of states 
and performing the integrals yields 


A = 


-E 

n^no 


{K \Hw\n){n\Hw\K°) 


niK — En 

g to) _ f 

"mx — En 


tb ta 


o-itf-ti)mK 


(4) 


- - ta)^K°\Hw\no){no\Hw\K°)e-^^f-^^^”^’< . 

The coefficient of 4 — ta is the desired Am ^. For states 
n with En < rriK there are exponentially growing con¬ 
tributions that complicate the extraction of the linear 
E — ta dependence. The control of these exponentially 
growing terms makes it clear that the methodology works 
best if the number of such contributions is small. A vari¬ 
ant of this general procedure can be used to compute 
long-distance contributions to 5 and rare Kaon decays 
mentioned in this work. Work along those lines is in 
progress [33) . 


THE RARE DECAY K ivuP 

The branching ratios for the and Kl ^ 

decays are given by 


BR(K+ 7r+z/iy(7)) = K+(l + Aem) 


A5 

A 


X{xt) 


Pe{X) + 




A5 


where 


-11 


= (5.173 ±0.025) • 10 

kl = (2.231 ±0.013) ■ 10“^^ 
Aem = -0.003 | 2 l] , 


A 


0.2252 

A 

0.2252 


(5) 

X{xt) ) (6) 

m , ( 7 ) 

m , ( 8 ) 

(9) 


X{xt) = 1.481 ± 0.005th ± 0.008. 
P.(A') = (0.404 ±0.024) 


'exp 
n 4 


, (10) 

[HIlIHlS] . (11) 


The theoretical results for these two decays read (keeping 
an explicit dependence on the CKM angles) [35] : 


BR{K+ TT+vv) = I (8.39 ± 0.30) x 10"“ 


iVbl 


n 2.8 


n 0.708 


73.2° 


40.7 X 10-3 

= I (3.36 ±0.05) x 


} , (12) 
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\Vub\ 1 

2 

iVbl 1 

2 

sin( 7 ) 17 

3.88 X 10-3 _ 


40.7 X 10-3 _ 


_sin(73.2o)J / 


}. (13) 


Note that the explicit CKM dependence in Eq. is 
accurate to about 1% in the ranges 37 x 10“3 < \Vch\ < 
45 X 10“3 and 60° < 7 < 80°. The non-CKM uncer¬ 
tainty on BR(Al ^ iT^uP) is much smaller than in the 
charged mode due to large uncertainties associated with 
the quantity Pc{X) introduced in Eq. ®. 

Using the complete unitarity triangle fit results to de¬ 
termine the relevant CKM entries, we obtain the follow¬ 
ing SM predictions: 


BR(A+ ^ TT+z/z/) = 

BR(Al ^ TT^Z/Z/) = 


f (8.64 ±0.60) X 10-^^ SM 
[(l7.3t}JJ) X 10-11 E949 [6] 

f (2.88 ±0.25) X 10-11 SM 
[<2.6x10-3 E391a [46] 


In the study presented below we consider two future sce¬ 
narios in which the experimental central value of the 
branching ratio either does not change 
or shifts to the SM prediction. The NA62 experiment 
should be able to reduce the uncertainties to (5exp = 7% 
and 10%, respectively. (These estimates are based on 
the expectation that NA62 will collect about 100 SM 
events by 2017 [H [9]; the 7% uncertainty is obtained 
by rescaling the expected number of events by the ratio 
BR(A+ ^ 7 r+z/z/)exp/BR(A+ ^ 7r+z/z/)sM.) The KOTO 
experiment at JPARC expects to observe ir^uP at 

the SM level with a 10% uncertainty 


DIRECT CP-VIOLATION IN AT ^ tttt AND e'/e 

The effective Hamiltonian responsible for contributions 
to e'je is, at scales larger than /i ^ 0{mc) [50] : 


y - 
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where Xq = Vq^V^^ and the current-current, penguin and 
semi-leptonic operators are: 


Ql = {SaC^)v-A{c^da)v-A , 

Q 2 = {sc)v-A{cd)v-A , 

Ql = {SaUl3)v-A{ui3da)v-A , 

Q 2 = {su)v-A{ud)v-A , 

(15) 

(16) 

(17) 

(18) 

Qs = {sd)v-A '^{qq)v-a , 

(19) 

Q 4 = {Sadp)v-A '^{QpQa)v-A , 

(20) 

Qs = {sd)v-A ^(^^)y+A , 

(21) 

Qe = {Sadp)v-A '^{QpQa)v+A , 

(22) 

Q 7 = ^{sd)v-A E ®g(99)v+A , 

(23) 

Qs = A(^ad/3)v-A y^,eq{qpqa)v^A 

, (24) 

Q 9 = 2 ^sd)v-A E *^q{d<l)v-A , 

(25) 

Qio = A(^adp)v-A y^,eq{qpqa)v-A 

, (26) 

where the sums in the operators Q 3-10 run over the four 
lightest quarks, denotes the quark electric charges, a 

and [3 denote color indices (omitted for color singlet op¬ 
erators), and {V ±A) = 7 ^( 1 ± 75 ). Below /i ^ 0{mc) the 
charm quark is integrated out and the effective Hamilto¬ 
nian is more commonly written as: 

Gf 

'HeS = E ± r Qi{y) 

i=l 

, (27) 


i 

(Qz)o 

{Qi)2 

1 

-0.151(44) 

0.00965(59) 

2 

0.169(56) 

0.00965(59) 

3 

-0.0492(661) 

0 

4 

0.271(111) 

0 

5 

-0.191(64) 

0 

6 

-0.379(128) 

0 

7 

0.219(61) 

0.286(11) 

8 

1.72(39) 

1.314(76) 

9 

-0.202(70) 

0.01447(89) 

10 

0.118(50) 

0.01447(89) 


TABLE L Current determinations of ^ (7r7r)/=o,2 matrix 
elements for each of the operators given in Eqs. (18 26). Eor 
7 = 0 and 7 = 2 we take /x = 1.531 GeV and // = 3 GeV, 
respectively. 


ryi = 1.87 ± 0.76 [59] mt,poie = (173.5 ± 1.0) GeV 
rys = 0.5765 ± 0.0065 [60] mc(mc) = (1.275 ± 0.025) GeV 
rj 3 = 0.496 ± 0.047 [SI] £exp = (2.232 ± 0.007) x 10“^ 
Bk = 0.766 ± 0.010 [IQ] /k = (156.3 ± 0.9) MeV [ID] 


TABLE 11. Inputs used in the calculation of £. The top and 
charm quark masses as well as quantities not explicitly given 
are taken from Refs. HUS]. 


where the quantities oj = ReA 2 /Re 74 o = 0.04454(12), 
Re^o = 3.3201(18) x 10“'^ GeV, Re^a = 1.4788(41) x 
10 -® GeV and 


where the operators Q1 2 have disappeared, the sums in 
the operators Q 3-10 run over the three lightest quarks, 
and r = -V/A„ = (note that 

An + Ac + At = 0 ). Complete NLO expressions for the 
coefficients Zi^ji) and ^i(/i) can be found in Refs. [5Ql[5T] . 

Following Refs. [5Ql [521455] . we write: 


Re 



Alo 

^2 


Re 


uj e 


^((52—(5o+7r/2) 




Im(A 2 ) Im(Ao) 
Re(A 2 ) Re(j4o) 


■ cos ((^2 - 5o + 7r/2 - (pe) 


V2\e\ 

'Im(^2) Im(Ao) 
Re(A 2 ) Re(j4o) 


(28) 


(29) 


Gp 


10 


El 2/i(M)] {{T^T^)l=0\Qi\K) , 

^ i=l 

(30) 


Gp 


10 


-^A„ E + T ViM] ((7r7r)/=2|Qi|-^), 

^ i=l 

(31) 


(pe = (43.5 ± 0.5)° (32) 

are taken from experiments HESj. The direct and ex¬ 
perimental determination of the phases 60^2 differ at the 
two sigma level 

^ =,5 _<5 1(42.3 ±1.5)° PDG 0 

^ ° 2 1(54-6 ±5.8)° RBG [571|58| ■ 

(33) 

Fortunately, due to the central value of the combination 
S 2 — So 7r/2 — (j)^ and to the large uncertainties in the 
determination of the various matrix elements, these two 
choices yield almost identical results; for definiteness, we 
follow the approach of Ref. [58] and use the phases ex¬ 
tracted from the lattice. The inclusion of isospin breaking 
corrections modifies the QCD penguins (Qs-e) contribu¬ 
tion to the (7r7r)/=o amplitude [62]: 

[Im(Ao)] 

QCDP ^ [Im(^o)]QCDP (4 ^eff) d ) (34) 

a = 1.017 ES] , (35) 

Ceff = (14.8 ± 8.0) X 10"2 IHUg . (36) 
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FIG. 1. Dependence of Re(£'/£) on fj. The horizontal yellow 
band is the la experimental measurement. The thin verti¬ 
cal red band is the Icr determination of fj from the standard 
unitarity triangle fit. The three lines correspond to taking 
ImAo = [ImAoJuBC + (0,1, 2) (5[Im^o]RBC- The shaded areas 
around these three lines are due to all the remaining sources 
of uncertainties (ImA 2 , fe, (5o, ^e). 


As we infer from Eq. (36), phenomenological estimates 
of isospin breaking effects vary over a considerable range 
and per Refs. [62H65] may be as high as about 15%. In 
any case these corrections are likely sub dominant to the 
uncertainties in the current lattice calculations [saisH] 
and are therefore being ignored at present. However, 
electromagnetic and isospin effects will be calculated in 
lattice QCD simulations in the next few years (See, for 
instance. Ref. [66] for a discussion of technical issues in¬ 
volved in these calculations). 

The imaginary part of the 1 = 0,2 matrix elements 
have been recently calculated and read 


Im(^2) = (-6.99 ± 0.20stat ± O.S^yat) x ^ 

(37) 

Im(Ao) = (-1.90 ± 1.23,tat ± l-OSsyst) x 10““ GeV . 

(38) 


The present experimental [STHzg and theoretical [58| 
results read: 


Re 1 

f ^exp 


\ ^exp 

Re 1 

(ik 

\ ^exp 



= (16.6 ± 2.3) X 10“^ , (39) 

= (1.36 ± 5.15atat ± 4.59syst) X 10-4 . (40) 


Where the notation clarifies that we normalize the the¬ 
oretical prediction for s' to the experimental measure¬ 
ment of 5 (that has a much smaller error than the cor¬ 
responding theory determination). For completeness we 
mention that the inclusion of the isospin breaking correc¬ 
tions given in Eqs. (|M|)-(36) shifts the SM prediction for 
Re (s'/s) to (0.5 ±5.9) x 10“^. The reduction in the total 




FIG. 2. Standard unitarity triangle (SUT) fit. In the top 
panel we show the impact of the e'/e measurement using the 
most recent calculations of the K (tttt) 7 = 0,2 matrix ele¬ 
ments. In the lower panel we consider a future scenario in 
which the uncertainty on ImAo and ImA 2 is reduced to 18% 
and 5%, respectively (see text for more details). The ImAo 
central value shifts to the value expected from the experimen¬ 
tal determination of s'/s and the standard unitarity triangle 
fit. All constraints are plotted at two two sigma level with the 
exception of ^exp/^exp for which we show both the one and 
two sigma contours. 


error is due to the factor a(l — f^eff) ^ 0.87 that multi¬ 
plies the QCD penguin contribution to ImAo, resulting 
in a difference of 2.5 a from the measured value. 

In passing we briefly note that the lattice QCD pre¬ 
diction in Eq. (|4Q| ) differs somewhat from the result pre¬ 
sented in Ref. [62 j . (1.9±4.5)xl0“^, which has a smaller 
total error. Consequently Ref. [62| gets a larger deviation 
of 2.9 a from the measured value of Re(5[j^/5exp) than 
the 2.1 cr that RBC-UKQCD [58| gets. This difference 
can be explained by a combination of (1) isospin break¬ 
ing corrections and (2) the use of the operator relation, 
Q 4 = Qs ± Q 2 — Qi, along with the experimental infor¬ 
mation on Re(Ao) and of an estimate of the ratio of the 
/ = 0 matrix elements of Qi and Q 2 (based on the RBC- 
UKQCD lattice QCD results and a large-calculation) 
to reduce the uncertainty on (Q 4 )o- Additionally, the au¬ 
thors of Ref. [62| used the RBC/UKQCD lattice result 
for the dominant {Q6)o contribution to Im(Ao)/Re(Ao), 
but expressed contributions proportional to (Q 3 , 5 , 7 - 10)0 
in terms of 1=0 matrix element ratios (by writing Re(Ao) 
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FIG. 3. Kaon unitarity triangle fits. Red contours (labelled KUT) are obtained including BR(iF^ ^ e^h/^exp, £ and 

\Vcb\- The small black contour is the current standard unitarity triangle (SUT) fit from B/K physics. In panel (a) we present 
the current status. The yellow area is allowed by sjh/^exp and the region below the blue curves is allowed by BR(i^^ ^ 7 r^z/z>). 
In panels (b-d) we show the impact of future improvements on the experimental determination of BR(iF^ ^ tt^z/z/) and on the 
theoretical calculation of the quantities ImAo and linA 2 (see text for more details). In panel (b) we assume that future central 
values for these quantities remain unchanged. In panel (c) we consider a scenario in which ImAo shifts to the value expected 
from the experimental determination of e'/e and the standard unitarity triangle fit. In panel (d) we assume, in addition, that 
the future experimental determination of BR(iF^ ^ tt^z/z/) will shift to the central value of the SM prediction. 


in terms of (Qi, 2 )o) that are then set to reference values 
inspired by large-A^. 

In both the / = 0 and 1 = 2 lattice computations, 
the systematic uncertainty is currently dominated by the 
perturbative truncation error in the computation of Wil¬ 
son coefficients and the matching from the RI to MS- 
bar scheme. This error is of the order where /r 

is the scale at which one matches to perturbation the¬ 
ory. By running non-perturbatively through the charm 
and eventually also bottom thresholds, these truncation 
errors will be reduced significantly in the near future. 
For /i = 50 GeV, e.g., perturbative truncation errors of 
0(1-2%) are feasible. Initial efforts along those lines are 
already in progress m- 

The ten / = 0 matrix elements are given in Table SII 


of Ref. [58] . The systematic uncertainties on the individ¬ 
ual matrix elements are obtained from those presented in 
Table II of Ref. [58| without including the “Wilson co¬ 
efficients” and “parametric errors” contributions and are 
about 22.5%. 

The 1 = 2 matrix elements in the continuum limit and 
in the ( 7 , 7 ) and (^,^) RI-SMOM schemes at 3 GeV are 
given in the last two rows of Table XIV of Ref. [57] . The 
conversion between the basis used in Ref. m and the 
standard operator basis we use is given in Eqs. (69-70) 
of Ref. [57] . The conversion of these matrix elements to 
the MS scheme at 3 GeV is achieved via the conversion 
matrix Cij = Sij -h with i^j = 1, 7,8. The 

relevant Avij entries are given in Table IX and XI of 
Ref. [ 75 ] for the (^,^) and ( 7 , 7 ) schemes, respectively. 
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FIG. 4. Impact of a future measurement of 

assuming SM central values with a 10% uncertainty [48|[49]. 


Uncertainties associated with RI-SMOM to MS conver¬ 
sion are estimated by comparing the results in the (^, ^) 
and ( 7 , 7 ) schemes, respectively. 

The matrix elements that we obtain are summarized in 
Table [I] (for / = 0 and / = 2 we take = 1.531 GeV and 
/i = 3 GeV, respectively), where, as discussed above, the 
systematic uncertainties on the / = 0 matrix elements 
are about 23% which have been added to the respective 
statistical uncertainties in quadrature. The uncertainties 
on the 1 = 2 matrix elements have been obtained by 
combining in quadrature the errors quoted in Table XIV 
of Ref. [57] with the scale uncertainty (the latter has 
been defined as the difference between the MS matrix 
elements obtained via an intermediate (^, or ( 7 , 7 ) RI- 
SMOM scheme). 

The remaining statistical and systematic uncertainties 
in the isospin symmetric limit are amenable to improve¬ 
ment by a straightforward numerical effort. We there¬ 
fore believe that an error on individual / = 0 operator 
matrix elements of order 5%-10% is achievable within 
five years given sufficient computational resources. Us¬ 
ing current central values of matrix elements, the corre¬ 
sponding propagated error on Im( 74 o) is 0 ( 10 - 20 %). 

Motivated by this discussion, we assume that the un¬ 
certainties on the / = 0 , 2 matrix elements will reduce to 
6 lm{A 2 ) = 5% and 6 lm{Ao) = 10%stat + 15%syst = 18% 
on a time-scale of five years. Note that these error es¬ 
timates are given with respect to the central values in 
Eqs. (37) and (38). Allowing for fluctuations of the cen¬ 
tral values of the matrix elements that control Im(Ao) 
the expected future theoretical errors on this quantity is 
^lm(Ao) 0.34 x 10 -^^ GeV. 

Before concluding this section let us comment on the 
treatment of the correlations between the errors on the 
matrix elements listed in Table [T| In the extraction of 


ImAj (/ = 0 , 2 ) given in Eqs. ( |37| ) and (38) all statisti¬ 
cal correlations have been taken into account; systematic 


uncertainties for ImAo were ascribed to individual oper¬ 
ator matrix elements and treated as uncorrelated (corre¬ 
lations tended to reduce the errors). On the other hand, 
we neglect correlations between the ImA 2 and ImAg er¬ 
rors. Given the difference in the uncertainties on these 
two matrix elements, correlations can be safely neglected 
at present but will need to be included when the total 
error on ImAo will become small enough. 


INDIRECT CP-VIOLATION AND s 

The basic expression for 5 is 

12^/27r2Am^P 
ViSoix,) + V2Soixt) 

+ 2v3So{Xe,Xt)VMVM , 


(41) 


where the numerical inputs we use are summarized in 
Table [TT| The quantity k,^ summarizes the impact of long 
distance effects and can be extracted from the knowl¬ 
edge of Im Ao and from an estimate of the long distance 
contributions to Am^- Eollowing Ref. m, we have: 

.. + ( 42 ) 

where p = 0.6 ± 0.3. Using the most recent RBC de¬ 


termination of Im(Ao) and (j)^ of Eq. (32), we obtain 
K,s = 0.963 ± 0.014 (See also the analysis presented in 
Ref. [77]). 

The single largest remaining parametric uncertainty on 
5 is due to I Vcb \ • The latter determines the parameter A 
of the CKM Wolfenstein parametrization that enters 5 
to the fourth power (via oc A^). Presently both 

inclusive and exclusive b ciu decays lead to a ^ 2% 
error on \ Vcb\] unfortunately a 2.5cr tension between these 
two determinations lead, via standard PDG rescaling, to 
a 2.5% error on the averaged result. Unless this tension 
signals some form of new physics, the most probable av¬ 
enue to improve the total uncertainty on 5 is hoping that 
future experimental and theoretical results for b —> ciu 
transitions will resolve this tension thus leading to a final 
sub-percent error on | Vcb \ • For completeness we mention 
that an alternative approach to unitarity triangle studies 
that does not make use of semileptonic decays has been 
proposed in Ref. [75] . 


RESULTS 

In Eig. we show the dependence of Re( 5 '/ 5 ) on 7 
for different choices of the Im(Ao) central value. Note 
how the uncertainty on this matrix element completely 
dominates the total uncertainty. 






















FIG. 5. Kaon unitarity triangle fits using e^h/^th rather than e^h/^exp- 


See the caption in Fig. for further details. 


In Fig. we present the standard unitarity trian¬ 
gle (SUT) fit obtained using B and K physics mea¬ 
surements. All the inputs used in the fit are taken 
from Refs. n nni Eu In particular, we adopt \Vcb\ = 
40.41(94) X 10“^ and \Vub\ = 3.68(26) x 10“^ obtained 
from a combination of inclusive and exclusive determi¬ 
nations of these CKM elements (because of the slight 
tensions between inclusive and exclusive {q = r^,c), 
the errors on the weighted averages are rescaled follow¬ 
ing the standard procedure described in Ref. ED- In 
the upper panel we show the present constraints imposed 
by (this notation underscores that in the theo¬ 

retical prediction for the e'/e ratio the denominator is 
taken from experiment and not calculated using the in¬ 
puts listed in Table 0 - In the lower panel we entertain 
a future scenario in which the uncertainty on ImAo is 
18%, as discussed above. Its central value is assumed to 
shift to what is necessary to reproduce the experimental 
determination of e'/e (note that even though e'/^exp is 
proportional to r/, the minimized over every other pa¬ 
rameter is not a symmetric function of fj). If the future 


ImAo central value does not shift, the e'/s allowed region 
is f] > 1.6 (see Fig. [^) implying a very strong tension 
with the standard fit. 

In Fig. we present the Kaon unitarity triangle fits 
(KUT) in various scenarios. In this fit we use only 
inputs from Kaon physics with the exception of tree- 
level determinations of \Vcb\ from inclusive and exclusive 
b c£u decays. The red contours are obtained including 
BR(iF+ ^ tt+z/F), e'/ s, 5, and \Vcb\ (from an average of 
inclusive and exclusive decays). The small black contour 
is the current standard unitarity triangle fit from B/K 
physics. 

In Fig.[^ we present the current status of this fit. The 
yellow area is allowed by e'/s and the region below the 
blue curves is allowed by BR(iF+ ^ tt+z^P). 

In Figs. [^-[^ we show the impact of future improve¬ 
ments on the experimental determination of BR(iF+ ^ 
tt+z^F) and on the theoretical calculation of the quan¬ 
tities ImAo and ImA 2 . In particular, in panel (b) we 
assume that future central values for these quantities re¬ 
main unchanged and that experimental and theoretical 
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FIG. 6. Impact of a future measurement of BR(i^L ^ 
assuming SM central values with a 10% uncertainty [48|[49]. 

uncertainties on BR(iF+ ^ and Im(74o,2) reduce 

as discussed above. In panel (c) we consider a scenario in 
which Im^do shifts to the value expected from ^exp/^exp 
and the standard unitarity triangle fit. In panel (d) we 
assume, in addition, that the future experimental deter¬ 
mination of BR(iF+ ^ TT+z/p) will shift to the central 
value of the SM prediction. 

Finally in Fig. we show the impact of a future mea¬ 
surement of BR{Kl at the SM level with 10% 

uncertainty. 

In Figs. and we present the very same fits but 
display rather than 5(.j^/5exp- The reason for con¬ 
sidering the theoretical ratio is that it is indepen¬ 

dent of the unique SM CP violating parameter fj (both 
numerator and denominator being proportional to f ]). As 
precise theoretical calculations of this ratio become avail¬ 
able, direct comparison with the corresponding experi¬ 
mental ratio can provide a new CP-conserving observable 
that can be used to constrain the relevant combination 
of CKM parameters in the UT-fit. 

CONCLUSIONS AND OUTLOOK 

In this work we have tried to draw attention to the sig¬ 
nificant progress that lattice QCD methods have recently 
made for quantitatively addressing non-perturbative ef¬ 
fects in several iF-decays such as iF ^ tttt and the direct 
CP-violation parameter Re {s'/e)^ the long-distance con¬ 
tribution to 5, and rare K-decays; of particular interest 
will be future lattice QCD determination of the long dis¬ 
tance contributions to the quantity Pc{X) that contains 
charm quark contributions to BR(iF+ ^ [80] . 

This means in the near future we will be able make 
better use of experimental data, existing and forthcom¬ 
ing, to better constrain the SM and search for new ef¬ 
fects. In particular, it appears that we can start to con¬ 


struct a unitarity triangle based primarily on K-physics. 
With improvements in the lattice calculations that are 
on the horizon and with results from forthcoming K- 
experiments a tighter KUT should soon become avail¬ 
able. It would be very valuable to compare the solution 
of such an improved KUT with the Standard Unitarity 
Triangle (SUT) coming primarily from 5-physics. 

In particular it now seems realistic that lattice calcu¬ 
lations can reduce the errors on Re(5Y5) to less than 
around 11% of the current experimental central value in 
about 5-years time. It may therefore be timely for the 
experimental community to plan an improved determi¬ 
nation of Re{e'/e)^ the current experimental errors on 
that quantity being around 15%. 


ACKNOWLEDGMENT 

C.L. and A.S. want to thank their collaborators from 
RBC-UKQCD for discussions and ecouragement. The 
work of CL and A.S. is supported in part by US DOE 
Contract #AC-02-98CH10886(BNL). The work of E.L. 
is supported in part by the United States Department of 
Energy under grant number DE-SC0010120. 


[1] N. Cabibbo, Meeting of the Italian Sehool of Physies 
and Weak Interactions Bologna, Italy, April 26-28, 1984, 
Phys. Rev. Lett., 10, 531 (1963), [,648(1963)]. 

[2] M. Kobayashi and T. Maskawa, Prog. Theor. Phys., 49, 
652 (1973) 

[3] J. Christenson, J. Cronin, V. Fitch, and R. Turlay, 
Phys.Rev.Lett., 13, 138 (1964) 

[4] K. Olive et al (Particle Data Group), Chin.Phys., C38, 
090001 (2014) 

[5] S. Adler et al. (E787), Phys.Rev.Lett., 79, 2204 (1997), 
arXiv:hep-ex/9708031 [hep-ex] 

[6] A. Artamonov et al. (E949), Phys.Rev.Lett., 101 , 191802 
(2008), arXiv:0808.2459 [hep-ex] 

[7] M. Moulson (NA62), (2013), arXiv:1310.7816 [hep-ex] 

[8] F. Newson, G. Aglieri Rinella, R. Aliberti, F. Ambrosino, 
B. Angelncci, et al, (2014), arXiv:1411.0109 [hep-ex] 

[9] A. Romano, (2014), arXiv: 1411.6546 [hep-ex] 

[10] S. Aoki, Y. Aoki, C. Bernard, T. Blum, G. Colangelo, 
et al., Eur.Phys.J., C74, 2890 (2014), arXiv:1310.8555 
[hep-lat] 

[11] A. J. Buras, M. E. Lautenbacher, and G. Ostermaier, 
Phys. Rev., D50, 3433 (1994), arXiv:hep-ph/9403384 
[hep-ph] 

[12] G. Buchalla and A. J. Buras, Phys. Lett., B333, 221 
(1994), arXiv:hep-ph/9405259 [hep-ph] 

[13] G. Buchalla and A. J. Buras, Phys. Rev., D54, 6782 
(1996), arXiv:hep-ph/9607447 [hep-ph] 

[14] A. J. Buras and R. Fleischer, Phys. Rev., D64, 115010 
(2001), arXiv:hep-ph/0104238 [hep-ph] 

[15] A. J. Buras, M. Gorbahn, U. Haisch, and U. Nierste, 
JHEP, 0611, 002 (2006), arXiv:hep-ph/0603079 [hep- 
ph] 




















10 


[16] E. Blucher, B. Winstein, and T. Yamanaka, 

[45] J. Brod and M. Gorbahn, Phys. Rev., D78, 034006 

Prog.Theor.Phys., 122, 81 (2009) 

(2008), arXiv:0805.4119 [hep-ph] 

[17] K. Agashe, G. Perez, and A. Soni, Phys.Rev., D71, 

[46] J. Ahn et al. (E391a), Phys.Rev., D81, 072004 (2010), 

016002 (2005), arXiv:hep-ph/0408134 [hep-ph] 

arXiv:0911.4789 [hep-ex] 

[18] M. Blanke, A. J. Buras, B. Duling, S. Gori, and 

[47] T. Yamanaka (KOTO), PTEP, 2012, 02B006 (2012) 

A. Weiler, JHEP, 0903, 001 (2009), arXiv:0809.1073 

[48] T. Komatsubara, Prog.Part.Nucl.Phys., 67, 995 (2012), 

[hep-ph] 

arXiv: 1203.6437 [hep-ex] 

[19] S. Casagrande, F. Goertz, U. Haisch, M. Neubert, and 

[49] K. Shiomi (KOTO), (2014), arXiv: 1411.4250 [hep-ex] 

T. Pfoh, JHEP, 1009, 014 (2010), arXiv: 1005.4315 [hep- 

[50] A. J. Buras, M. Jamin, and M. E. Lautenbacher, 

ph] 

Nucl.Phys., B408, 209 (1993), arXiv:hep-ph/9303284 

[20] L. S. Littenberg, Phys.Rev., D39, 3322 (1989) 

[21] J. A. Bailey et al. (MILG), Phys. Rev., D92, 034506 

[hep-ph] 

[51] G. Buchalla, A. J. Buras, and M. E. Lautenbacher, 

(2015), arXiv: 1503.07237 [hep-lat] 

Rev.Mod.Phys., 68, 1125 (1996), arXiv:hep-ph/9512380 

[22] X. Eeng, N. H. Ghrist, A. Portelli, and G. Sachrajda, 
PoS, LATTICE2014, 367 (2015). 

[hep-ph] 

[52] M. Giuchini, E. Eranco, G. Martinelli, L. Reina, and 

[23] A. J. Buras, M. E. Lautenbacher, M. Misiak, and 

L. Silvestrini, Z. Phys., C68, 239 (1995), arXiv:hep- 

M. Munz, Nucl.Phys., B423, 349 (1994), arXiv:hep- 

ph/9501265 [hep-ph] 

ph/9402347 [hep-ph] 

[53] A. J. Buras and L. Silvestrini, Nucl.Phys., B569, 3 

[24] H. Greenlee, Phys.Rev., D42, 3724 (1990) 

[25] We want to thank Laurence Littenberg for discussions 

(2000), arXiv:hep-ph/9812392 [hep-ph] 

[54] S. Bosch, A. Buras, M. Gorbahn, S. Jager, M. Jamin, 

regarding this background. 

[26] G. Beall, M. Bander, and A. Soni, 11th International 

et al, Nucl.Phys., B565, 3 (2000), arXiv:hep-ph/9904408 
[hep-ph] 

Symposium on Lepton and Photon Interaetions at High 

[55] A. J. Buras and M. Jamin, JHEP, 0401, 048 (2004), 

Energies Ithaea, New York, August 1-9, 1983, Phys. Rev. 
Lett., 48, 848 (1982) 

arXiv:hep-ph/0306217 [hep-ph] 

[56] Q. Liu, Ph. D. thesis, Golumbia University (2012). 

[27] L. Lellouch and M. Luscher, Gommun.Math.Phys., 219, 

[57] T. Blum, P. Boyle, N. Ghrist, J. Prison, N. Garron, et al, 

31 (2001), arXiv:hep-lat/0003023 [hep-lat] 

Phys.Rev., D91, 074502 (2015), arXiv:1502.00263 [hep- 

[28] G. Kim, G. Sachrajda, and S. R. Sharpe, Nucl.Phys., 

lat] 

B727, 218 (2005), arXiv:hep-lat/0507006 [hep-lat] 

[58] Z. Bai, T. Blum, P. Boyle, N. Ghrist, J. Prison, et al., 

[29] R. A. Briceo, PoS, LATTICE2014, 008 (2015), 

(2015), arXiv:1505.07863 [hep-lat] 

arXiv: 1411.6944 [hep-lat] 

[59] J. Brod and M. Gorbahn, Phys.Rev.Lett., 108, 121801 

[30] N. H. Ghrist (RBG, UKQGD), (2010), arXiv:1012.6034 

(2012), arXiv:1108.2036 [hep-ph] 

[hep-lat], 

[60] A. J. Buras, M. Jamin, and P. H. Weisz, Nucl.Phys., 

[31] Z. Bai, N. Ghrist, T. Izubuchi, G. Sachrajda, 

B347, 491 (1990) 

A. Soni, et al, Phys.Rev.Lett., 113, 112003 (2014), 

[61] J. Brod and M. Gorbahn, Phys.Rev., D82, 094026 

arXiv:1406.0916 [hep-lat] 

[32] N. H. Ghrist, X. Eeng, G. Martinelli, and G. T. Sachra¬ 

(2010), arXiv: 1007.0684 [hep-ph] 

[62] A. J. Buras, M. Gorbahn, S. Jager, and M. Jamin, 

jda, Phys.Rev., D91, 114510 (2015), arXiv:1504.01170 

(2015), arXiv:1507.06345 [hep-ph] 

[hep-lat], 

[63] V. Girigliano, G. Ecker, H. Neufeld, and A. Pich, 

[33] N. Ghrist, T. Izubuchi, G. T. Sachrajda, A. Soni, 

Eur.Phys.J., C33, 369 (2004), arXiv:hep-ph/0310351 

and J. Yu (RBG, UKQGD), PoS, LATTICE2013, 397 

[hep-ph] 

(2014), arXiv: 1402.2577 [hep-lat] 

[64] A. Buras and J. Gerard, Phys.Lett., B192, 156 (1987) 

[34] U. Haisch, PoS, KAON, 056 (2008), arXiv:0707.3098 

[65] V. Girigliano, G. Ecker, H. Neufeld, and A. Pich, 

[hep-ph] 

Phys. Rev. Lett., 91, 162001 (2003), arXiv:hep- 

[35] A. J. Buras, D. Buttazzo, J. Girrbach-Noe, and R. Kneg- 

ph/0307030 [hep-ph] 

jens, (2015), arXiv: 1503.02693 [hep-ph] 

[36] R. Knegjens, (2015), arXiv:1505.04928 [hep-ph] 

[66] N. Garrasco, V. Lubicz, G. Martinelli, G. T. Sachrajda, 
N. Tantalo, G. Tarantino, and M. Testa, Phys. Rev., 

[37] E. Mescia and G. Smith, Phys.Rev., D76, 034017 (2007), 

D91, 074506 (2015), arXiv: 1502.00257 [hep-lat] 

arXiv:0705.2025 [hep-ph] 

[67] V. Eanti et al. (NA48), Phys.Lett., B465, 335 (1999), 

[38] G. Buchalla and A. J. Buras, Nucl.Phys., B400, 225 

arXiv:hep-ex/9909022 [hep-ex] 

(1993) 

[68] A. Lai et al. (NA48), Eur.Phys.J., C22, 231 (2001), 

[39] M. Misiak and J. Urban, Phys.Lett., B451, 161 (1999), 

arXiv:hep-ex/0110019 [hep-ex] 

arXiv:hep-ph/9901278 [hep-ph] 

[69] J. Batley et al. (NA48), Phys.Lett., B544, 97 (2002), 

[40] J. Brod, M. Gorbahn, and E. Stamou, Phys.Rev., D83, 

arXiv:hep-ex/0208009 [hep-ex] 

034030 (2011), arXiv: 1009.0947 [hep-ph] 

[70] A. Alavi-Harati et al. (KTeV), Phys.Rev.Lett., 83, 22 

[41] G. Buchalla and A. J. Buras, Nucl.Phys., B548, 309 

(1999), arXiv:hep-ex/9905060 [hep-ex] 

(1999), arXiv:hep-ph/9901288 [hep-ph] 

[71] A. Alavi-Harati et al. (KTeV), Phys.Rev., D67, 012005 

[42] G. Buchalla and A. J. Buras, Nucl.Phys., B412, 106 

(2003), arXiv:hep-ex/0208007 [hep-ex] 

(1994), arXiv:hep-ph/9308272 [hep-ph] 

[43] A. Buras, M. Gorbahn, U. Haisch, and U. Nier- 

[72] G. Barr et al. (NA31), Phys.Lett., B317, 233 (1993) 

[73] L. Gibbons, A. Barker, R. A. Briere, G. Makoff, V. Pa- 

ste, Phys. Rev. Lett., 95, 261805 (2005), arXiv:hep- 

padimitriou, et al., Phys.Rev.Lett., 70, 1203 (1993). 

ph/0508165 [hep-ph] 

[74] J. Prison, P. Boyle, and N. Garron, PoS, LAT- 

[44] G. Isidori, E. Mescia, and G. Smith, Nucl.Phys., B718, 

TICE2014, 285 (2015), arXiv: 1412.0834 [hep-lat] 

319 (2005), arXiv:hep-ph/0503107 [hep-ph] 









11 


[75] C. Lehner and C. Sturm, Phys.Rev., D84, 014001 (2011), 
arXiv: 1104.4948 [hep-ph] 

[76] A. J. Buras, D. Guadagnoli, and G. Isidori, Phys.Lett., 
B688, 309 (2010), arXiv:1002.3612 [hep-ph] 

[77] Z. Ligeti and F. Sala, (2016), arXiv: 1602.08494 [hep-ph] 


[78] E. Lunghi and A. Soni, Phys. Rev. Lett., 104, 251802 
(2010), arXiv:0912.0002 [hep-ph] 

[79] Y. Amhis et al. (Heavy Flavor Averaging Group 
(HFAG)), (2014), arXiv:1412.7515 [hep-ex] 

[80] G. Isidori, G. Martinelli, and P. Turchetti, Phys. Lett., 
B633, 75 (2006), arXiv:hep-lat/0506026 [hep-lat] 


